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t>- ' Abstract 

(N ; 

In Refs. [21 E] two different approaches have been developed for the numerical 
computation of the effective energy in the presence of a magnetic flux tube, by using 
qq ■ the phase shift method. However, the opinion has been expressed that these two 

variations of the phase shift method are not equivalent. In this paper we aim to solve 
this ambiguity by comparing these two different approaches and showing that they give 
identical numerical results, within the numerical accuracy of the method. 
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^ ' 1 Introduction 
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The numerical computation of the fermion- induced effective energy 1 , in the presence of 
inhomogeneous magnetic fields of the form of a flux tube, is a topic that has attract the 
>. attention of several authors, see Refs. [TJ E3 E| . 

In Refs. E] two different approaches have been developed for the numerical computa- 
tion of the quantum energy in the presence of a magnetic flux tube, by using the phase shift 
method. The question that arises is whether these two different approaches are equivalent or 
not. In this paper we show, mainly by performing numerical computations, the equivalence 
of these two different approaches. In this way we give a detailed answer to the objections 
in Ref. [I], according to which the integral 2 we used for the numerical computation of the 
3+1 effective energy is not properly renormalized, and as a result it does not give the correct 
values for the effective energy. 

For this reason, in section 2 we present in a detailed way the renormalization procedure 
we followed in Ref. [2], and we emphasize that exactly the same method has been used by 
other authors, and especially in Ref. jHj. In section 3 we present our numerical results for 
the integrand of Eq. (12) for the quantum energy, in the case of the Gaussian magnetic 
flux tube (of Eq. (16) below). These results indicate a highly convergent behavior for the 
integral of Eq. (12). In section 4, we compare our results (Fig. El in this work), for the 
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lr rhe effective energy can be written as -E e // = E c i as + Eq, where E c i as is the classical energy and Eq is 
the quantum part of the fermion induced effective energy. 
2 of Eq. (26) in Ref. or in Eq. (12) in this work 
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quantum energy in the presence of the Gaussian flux tube, with the correspondent results 
of Ref. JTj (Fig. 4 in Ref. [4 ), and we obtain that they are identical, within the numerical 
error of the method. Finally, in section 5 we compare the two methods more carefully by 
using a different way (comparing with the vacuum polarization diagram) and we find again 
that they should be equivalent. 



2 Renormalization of the quantum energy 

In this section we explain in detail the renormalization procedure we used in Ref. 0. 

The quantum part of the fermion-induced effective energy in the presence of a magnetic 
field is given by the equation: 

E Q = -^Trln(^ + m}) (1) 

where fl = 7^-D M (/x = 0, 1,2,3) and = <9 M — ieA^. The gamma matrices satisfy the 
relationship {7^,7^} = 2g M1/ , and T is the total length of time. 

We have shown in Ref. |2 a that the 2 + 1 quantum energy for a magnetic field of the form 
of a flux tube is given by the equation 

1 r+°° k 

£q( 2+ i) = 5- / dk (A(fc) - c) (2) 

2vr Jo Jk 2 + m ) 

where c = limfc^ +00 A(k). Our numerical study, in Ref. |2j, shows that c = — 7r0 2 (0 = 
e$/27r). This means that c is independent of the special form of the magnetic field we 
examine and depends only on the total magnetic flux of the field $. We also note that c 
arises from an integration by parts (for details see Ref. 0) and it is not a quantity we 
subtract and add back in order to make the above integral convergent. 

The function A(k) was defined in Ref. j2] as the symmetrical sum over 1: 

A(*)= lim £ *,,.(*) (3) 

s,l=—L 

where Si tS (k) is the phase shift which corresponds to I th partial wave with momentum k and 
spin s. 

It is well known that there are no divergencies in 2+1 dimensional QED. From this 
point of view we expect the above integral of Eq. (2) to be convergent. In addition, our 
numerical results, for the integrand of Eq. (2), confirm with very good accuracy this expected 
convergent behavior, see Fig. 1 in Ref. [2] , or Fig. and Fig. El below in this paper. 
The calculation of the phase shifts have been performed by solving an ordinary differential 
equation. For details see Refs. [21 Ej. 

The corresponding 3 + 1 dimensional result is obtained by the 2 + 1 quantum energy if 
we replace mj by k% + rnj in (2) and integrate over the k$ momentum j^]. An overall factor 
of 2 (from the Dirac trace) must also be included. 

1 /-+oo /-A/2 flh n U 

e qm =- dh L, o it /,„ o( a w - c ) ( 4 ) 



n Jo ' J -a/2 2n/L z ^k 2 + kj + 
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As the integral over k$ is divergent we have introduced a regularization parameter A/2. Now, 
if we integrate over we obtain the unrenormalized quantum energy per unit length L z (L z 
is the length of the space box towards z direction) equal to 

1 r+°° I ' k 2 + m 2 f\ 1 A 2 f+°° 

In the weak field limit eB « m 2 , or for large rrif, the above expression for the effective 
energy must tend to the first diagram (vacuum polarization diagram) of the perturbative 
expansion of the effective energy: 



y Q(3+i) - g7r 3 j "^v- y o ^yyi-vyyi m 2 



+ -J— ln(^l) / d 2 xB 2 (x) 
24vr 2 m 2; i V ; 



(6) 



where B(q) = J d 2 x e~ l ^' x B(x), and we have assumed the rescaling B — > e£?. 

Note that the vacuum polarization diagram for the 3+1 dimensional case (see Eq. (6) 
above) can be obtained from 2+1 dimensional result 

«( 2+1 ) 16 W yi Wl 7o Jm 2 +x(l-x)q 2 K) 



if we set m 2 — > /c 2 + m 2 and integrate over the &3 momentum (an overall factor of 2 must 
also be included). 

By comparing Eqs. (5) and (6) in the weak field limit we obtain that 

r+oo 1 r 

/ dkk(A(k) - c ) = — d 2 xB 2 (x) (8) 
Jo 12 J 

Note that we have confirmed the above equation numerically by performing numerical com- 
putation for several magnetic fields of the form of a flux tube. 

Now Eq. (5), for the unrenormalized quantum energy, can be written as 

2 i ™2 1 



1 f +0 ° (k 2 + mi\ , , s 



247r 2 m 



The logarithmically divergent term can be incorporated in the classical energy, as is shown 
below 



Eeff — Eclass + E, 

in i , , a 2 

2 W 



i 

1 r+oo ( k 2 + m 2 f \ 

+ 2^-1 ^)(A W -e) (10) 
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In this way the logarithmically divergent part corresponds to a charge renormalization ac- 
cording to the equation 

1 1 1 , / A \ 
and thus the renormalized quantum energy is 

E< ^ = -^sr dmn ^f) [Mk) - c) <i2) 

The above presented renormalization procedure corresponds to the on-shell renormalization 
condition 11(0) = (for more details see Ref. [Sj). 

Now, as the divergent part has been removed, the integral of Eq. (12), for the renormalized 
quantum energy, is expected to be convergent. In the next section, we will indeed see that 
our numerical results indicate a highly convergent behavior for this integral. 

Note that exactly the same renormalization procedure has been used by the authors of 
Ref. Particularly, in Ref. |Hj the analytical result for the 2+1 quantum energy for a 
special form of an inhomogeneous magnetic field is taken for granted from their previous work 
of Ref PU]- Then the replacement m'j — > k\ + m'j is performed, and the 3+1 dimensional 
result is obtained by integrating out the k% momentum. As we did in this paper, a cut- 
off A/2 is used in Ref. [Hj for the regularization of the integral over k%, and the same 
logarithmically divergent term arises. This term can be removed as it corresponds to a 
charge renormalization. In Ref. [S] we see that the remaining part, which corresponds to 
the regularized quantum energy, is an analytical formula with no divergencies. 



3 Convergence of the integral of Eq. (12) 

For the numerical computations we have chosen the Gaussian magnetic field 

B a (r) = ^ e*p(~) (13) 

where d is the spatial size of the magnetic flux tube, <fi = e$/2-7r, and $ is the total magnetic 
flux of the field. 

In Fig. [T]we have plotted the integrand of Eq. (12) 

I(k) = 2^1* (^p) ( A (*) - c ) ( 14 ) 

as a function of k for = 2.5, 4.5, 6.5, 8.5, d = 1 and rrif = 1. 

We see that our numerical results in Fig. indicate a highly convergent behavior for the 
integral of Eq. (12), as the integrand I(k) becomes zero for sufficiently large values of k. 



4 Comparing with the results of Ref. |3| 

As the convergence of the integral of Eq. (12) is not sufficient to show that our numerical 
results for the quantum energy are the same with those of the phase shift approach of Ref. [3] , 
we should perform a comparison of the numerical results of these two different approaches. 
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Figure 1: The integrand I(k) as a function of k for cf) = 2.5, 4.5, 6.5, 8.5 and m/ = 1, in the 
case of the Gaussian magnetic field of Eq. (16). 

For this we have plotted Eq^^^/ <p 2 as a function of d for <ft = 2.5, 4.5, 6.5, 8.5 and mj = 1, 
in the right-hand panel of Fig. El We see that our results in this figure are identical (within 
the numerical accuracy of the method) with the corresponding results of Ref. [I] (see the 
right-hand panel of Fig. (4) in Ref. 

Also, it is worth to compare the results of the phase shift method with the derivative 
expansion, which is the standard approximative tool in the case of smooth inhomogeneous 
magnetic fields (see for example Ref. [S]). In Fig. E]we see that our results are in a very good 
agreement with the corresponding results of the derivative expansion in the homogeneous 
limit \j \JB m << d, as is expected, where B m = 20/ 'd 2 . 
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Figure 2: Derivative Expansion versus Phase Shift Method for the Gaussian magnetic flux 
tube, for = 2.5, 4.5, 6.5, 8.5 and my = 1. 

5 Phase shift approach and vacuum polarization dia- 
gram 

In Ref. j3] the renormalized quantum energy 3 Eq, z+X x is computed as a sum of two terms 

E Q(3+1) = E ph(3+1) + E Q(3+1) ( 15 ) 

where the first term corresponds to the contribution of the phase shifts 

1 r+°° ( k 2 + m 2 f \ 

E ph{3+1) = -— 2 j o dkkln [-^T 1 ) (16) 

(2) 

and the second term Eq[^ +1 ^ is the renormalized quantum energy which corresponds to the 
vacuum polarization diagram (see Eq. (8)). 
The function A(k) is defined as 

£(*)= fim £ \s{k) (17) 

s,l=—L 

where 

1jM = 5iM - (k) - <A (2 j (*) (18) 

3 In what will follow we have dropped the index (ren) from the renormalized quantum energy Eq?£Ii\- 




- A(k)-c 

■ A(kH 2 A (2) (k) 
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Figure 3: A(k) — c and A(A;) = A(k) — cf) 2 A^ 2 \k) as a function of k for = 6.5, mj = 1 and 
d=l, where c = — 7r0 2 . 

and 8\g (k), Sl 2 ](k) are the coefficients of and 2 which arise if we expand the phase shift 
Si tS (k) in powers of 0. 

We emphasize that d^ 1 ] (k) and 5^ 2 ](k) are not the first and second Born approximations, 
as the Born expansion is an expansion in powers of the potential v^ s (r) 4 of the corresponding 
scattering problem, and this potential includes both terms of and 2 . 

From Eqs. (17) and (18) we obtain 

~A(k) = A(k) - 2 A (2) (A;) (19) 

as we can prove, from the integral formula for the first Born approximation, that A^^A;) = 0. 
Note that the functions A^(k) and A^(fc) are defined as: 



L L 




In Fig. EJwe have plotted A(k) - c = A(k) + tt0 2 and A(fc) = A(k) - <p 2 A {2) (k) as a 
function of k for = 6.5, rrif = 1 and d=l. We see that very shortly (k > 2) these two 
functions become identical. The main reason is that the function A^^(fc) tents rapidly to 
—7i, as we see in the left-hand panel of Fig. 0] 

4 For the exact form of the potential in the presence of the flux tube see Eq. (17) in Ref. [2] 
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Figure 4: A^^k) + n and i2(^) as a function of k for d=l and rrif 
A^ 2 \k) -> -7r as fc -> +oo. 



1. We see that 



The immediate consequence of the above discussion is that the integrals of Eqs. (12) 
and (16) have exactly the same convergent properties. Thus, Eq. (16), which was used for 
the computation of the quantum energy in Ref . [3] , has no additional advantage against Eq. 
(12) which we used in our work of Ref. [2]. On the other hand, if we use Eq. (12) for the 
quantum energy, it is not necessary to compute additional quantities like the phase shifts 
&f s }(k) and the vacuum polarization diagram. 

In addition, in the previous section we see that the results of these two variations of the 
phase shift method are in a very good agreement, however in this section we will go over 
this topic more carefully. 

If we demand the results of two methods to be equal (Eq@+i) = E'q, 3+1 \ see Eqs (12) 
and (15)), from Eqs. (6), (12), (15), (16) we obtain 



2tt 2 
1 



+°° (k 2 + mV 
dkk In k — - 



o 



m 



(A(k) 



37T 3 JO 



dxx(l 



X 



f 

oo 



2tt 2 



dkk In 



k 2 + m 2 f \ 



m. 



A(k) 



B{q)\ 2 ln{^- 



+ q 2 x(l — x) 



m, 



(20) 
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Figure 5: We have plotted the first and the second member of Eq. (21), multiplied by 
—d 2 /<f) 2 , as a function of d for rrif = 1. The continuous line corresponds to the vacuum 
polarization diagram, and the discrete points to the phase shift integral. According to this 
figure, Eq. (21) is valid. 

or, if we use Eq. (19), we obtain the equivalent equation 

2tt z Jo \ mj J 

1 r 1 r+°° ~ „ m 2 f + q 2 x(l — x) 

= -— dxx(l-x) dqq\B(q)\ 2 \n( f ^ >-) (21) 

07T 3 Jo Jo mj 

In Fig. Elwe have plotted the first and the second member of Eq. (21), multiplied by — d 2 /(f) 2 , 
as a function of d for m/ = 1. The continuous line corresponds to the vacuum polarization 
diagram, and the discrete points to the phase shift integral of the first member of Eq. (21). 
Thus, from Fig. it is seen that the first and second member of Eq. (21) are equal (the 
deviations between the numerical values of the first and second member of Eq. (21) are 
of the order of 0.1 per cent). As a result the two alternative approaches, for the quantum 
energy in the presence of a flux tube, are equivalent. 

In addition, if take into account the above results, we can compute the difference \Eqi^ + \\ — 
Eq, 3+1 \\, where Eq(s+i) and E'q^ +1 ^ are the renormalized quantum energies which correspond 
to the two variations of the phase shift method, see also Eq. (12) and (16). We see that 
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\Eq(3+i) — E'q(3+i)\/ Eq(3+i) ~ 0.001, which is of the order of the accuracy of the numerical 
computations, and as a result the numerical values for Eq^ + i) and Eq^ 3+1 ^ are identical. 

Note that, according to the above discussion, the first member of Eq. (21) can be viewed 
as the phase shift representation of the vacuum polarization diagram, or we can write 

= -^C dkkla f^r) < 22 > 

In order to study the convergent properties of the above integral, we have plotted, in the 
right-hand panel of Fig. the integrand of Eq. (22) 

«*>=2^ k p^r)( A<2) «H (23) 

for d = 1 and m/ = 1. We see that that Fig. 0] indicates a highly convergent behavior for 
the integral of Eq. (22). 

6 Conclusions 

In this paper we compared the results of the two different variations of the phase shift method, 
of Refs. E] , for the numerical computation of the quantum energy in the presence of a 
magnetic flux tube. In section 3 and 4 we show, mainly by performing numerical compu- 
tations, that these two alternative approaches give identical values for the quantum energy 
(|-Eq(3 + i) — Eq, 3+1 \\/Eq( 3+ i) ~ 0.001, see section 4), and as a result they are equivalent. 

In addition, in section 4 we show that in order to make the integral over k, of Eq. (12), 
convergent, it is not necessary to subtract the and 2 parts of the phase shift 5[ tS (k) (see 
Eq. 18) and to add them back in their Feynman diagrammatic form, as is done in Ref. [3]. 
This immediate convergence of the integral of Eq. (12) is a consequence of the translation 
invariance along the £3 axis. However, in the case of problems, e.g. with spherical symmetry, 
where the translation invariance is violated, the subtraction of the asymptotic part of the 
phase shift sum Y^ti^o $i,s{k) ( or t ne subtraction of the first and second Born approximation), 
in order to make the integral over k convergent, is unavoidable (see Refs. [§1 ITT]). 

It is worth to mention, that for the computation of the quantum energy in the presence 
of a magnetic flux tube, two other alternative methods have been developed: a) the Jost 
function method in Ref. pQ and b) the method of the worldline numerics in Ref. [Sj. Note, 
that in Ref. |2] we have compared our results to the Jost function method, in the case of a 
discontinuous magnetic field which is constant inside a cylinder of radius d, and zero outside 
it. From Fig. 4 in Ref. j2] and the Fig. 3 in Ref. pQ we see a very good agreement between 
the phase shift and Jost function method. 

Finally, it would be interesting if the two alternative approaches ( the Jost function and 
the worldline method) could give a figure like that of the right-hand panel of Fig. El of this 
work, for a smooth realistic magnetic field like the Gaussian magnetic flux tube. 
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